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In a class of multidimensional models, topology of the thick brane provides 
three chiral fermionic families with hierarchical masses and mixings in the ef- 
fective four-dimensional theory, while the full model contains a single vector- 
like generation. We discuss how to incorporate three non-degenerate neutrino 
masses in these models with the help of only one singlet bulk fermion. 
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1 Introduction 



One of the interesting possibilities which open up in theories with more than four 
spacetime dimensions is to explain the misterious pattern of fermion mass hierarchies 
© S 0] • m previous works [fj, |3|, we have shown how a single family of fermions, 
with vectorlike couplings to the Standard Model group 577(3) x SU(2) x U{1) in 
6 dimensions was reduced to n chiral families in 4 dimensions. This mechanism 
hinged on a localisation on a vortex in the "transverse" dimensions, characterized by 
a winding number n. As a variant, an effective vortex achieving the same result could 
be simulated by coupling the fermions to a winding-number 1 scalar elevated to the 
nth power. 

The value n = 3 is not favoured in such schemes, and was only invoked for its 
phenomenological interest. We also showed how a relatively simple scalar structure 
(the usual Higgs doublet H, supplemented by the vortex-defining field $) suffices to 
generate the mass hierarchy of the n families, while an auxiliary field X could be used 
to generate a mixing taking place mostly between adjacent generations. 

The key to the construction is that interaction with a vortex with winding number 
n leads, by the index theorem || , to n chiral zero modes in four dimensions, and this 
for each fermion species coupled to the structure. 

We now return to this construction and pay special attention to the case of neu- 
trinos. The simplest attitude (which was in a way implicit in the previous papers) 
would be to treat neutrinos just like the other fermions, that is to include both left- 
and right-handed fields, say v and N in four dimensions, and to generalize them to six 
dimensions. Dirac masses and mixings are then obtained just like those of the charged 
fermions. There is nothing fundamental against this approach, except the very small 
value of the ratio between neutrino and charged fermion masses (typically less than 
10~ 6 ). This would require severe fine tuning, either directly at the level of the Yukawa 
couplings themselves, or in a more covert fashion through a more complicated scalar 
structure (see footnote below). 

We therefore try to explore other ways of getting small neutrino masses. 

It has become standard practice || in the context of extra (flat and compacti- 
fied) dimensions, to assume that right-handed neutrions, being gauge singlets, could 
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propagate in the "bulk" of space, escaping confinement^. This is of course possible 
since such particles do not participate directly in the gauge interactions, and there- 
fore do not affect (at least at tree level, see for instance Ref. 0) the behaviour of 
electroweak exchanges. In such a context, the smallness of the neutrino mass stems 
not from artficially tuned Yukawa couplings, but from the severely reduced overlap 
between the confined and "bulk" wave functions. Could such a strategy work here? 
The transcription of such a mechanism in our context involves the introduction of a 
six-dimensional fermion field, N, uncoupled to the vortex field $. This however is 
in general not sufficient, if the two-dimensional transverse space stays unbounded (a 
point which we did not really need to specify this far, since the vortex achieved in 
any case the localisation of physical fields in a four- dimensional tube). Using 9 and 
r as variables to describe the two extra dimensions, we should now require r < R, 
while keeping R 3> r vortex , to avoid perturbing the previous construction. Is this 
step sufficient? In principle, if 1/R were close to the expected light neutrino mass 
scale (say a fraction of an eV), we would have a whole set of Kaluza-Klein towers at 
hand, each with relatively light states to couple to the n left-handed neutrinos. This 
line is however somewhat dangerous, since it may lead to a conflict with data on the 
evolution of supernovae f|. We will avoid this situation by requesting 1/R ^ m u . 
For simplicity we may even take[] 1/R^> 100 GeV. 

This further requirement however leads us to a pecular situation: we may not 
have enough light "right-handed" partners to provide masses to the n neutrinos! 
Indeed, we typically expect only one fundamental mode of zero mass. Introducing 
n such fermions would of course solve the problem, but would be in opposition to 
the approach, where at least the charged family multiplication simply results from 
the topological structure. Fortunately, even with the simplest use of one N field, we 
get easily two light states (this is a funciton of the chosen charge assignements, see 

* Other possibilities include putting left- and right-handed neutrinos to very different places inside 
the thick brane ||] or putting the scalar (not the fermion) in the bulk. The latter approach can be 
easily implemented in the model of Ref. || where additional scalar fields, others than the electroweak 
Higgs doublet, participate in mass generation. 

^This is to avoid problems with weak interactions universality. We suppose that the Standard 
Model gauge fields are localized inside the brane by means of a mechanism [^|, |l^| which preserves 
charge universality by itself, see discussion in Ref. Jll[] . 
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below), leaving only n — 2 degenerate massless vt states. Since one strictly massless 
neutrino is no problem (only differences in masses squared are tested this far), the 
neutrino degeneracy is thus solved for n < 3. 

The scheme presented below uses one extra scalar field to match quantum numbers 
and generate the neutrino masses. It is of course also possible to induce mixing terms, 
but these are probably redundant, since off-diagonal terms are already allowed for in 
the charged leptons sector. 



2 Models with a single generation in extra dimen- 



If different fermionic modes have different wave function profiles in extra dimensions, 
then their overlaps with the Higgs wave function may produce hierarchical structure 
of masses and mixings In the class of models [0, each multi-dimensional 
fermion develops three chiral zero modes localized on a four-dimensional brane due 
to topological properties of the brane background. The Index theorem guarantees 
that the three zero modes are linearly independent, and thus have different profiles in 
extra dimensions. Analysis of the equations for these zero modes demonstrates that 
a hierarchy in the mass matrix indeed appears due to overlaps of the wave functions. 
For the discussion of this mechanism and comparison with other approaches, see 



realisation of this approach which was suggested in Ref. ||. The model is formulated 
in six dimensions, and the brane is an Abelian vortex made of the gauge field of U(l) g 
gauge group and a scalar field Field content of the model of Ref. [||] is given in 
Table 1 for easy reference (we concentrate here on the leptonic sector since we are 
interested in neutrino masses). Like other Standard Model fermions, leptons are zero 
modes of the six-dimensional Dirac spinors, namely, of the electroweak SU(2) doublet 
L and SU(2) singlet E which have axial charges +3/2 and —3/2, respectively, under 
U(l) g . The interaction with vortex field, 
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fields 




profiles 


charges 


representations 










U(l)y 


SU(2) W 


SU(3) C 


scalar 


$ 


F(r)e ie 
F(0) = 0, F{oo) = Vq, 


+1 





1 


1 


scalar 


X 


X{r) 

X(0) = v x , X(oo) = 


+1 





1 


1 


scalar 


H 


H(r) 

H(0) = v H , H(oo) = 


-1 


+ 1/2 


2 


1 


fermion 


L 


3 L zero modes 


axial +3/2 


-1/2 


2 


1 


fermion 


E 


3 R zero modes 


axial —3/2 


-1 


1 


1 



Table 1: Field content of the model of Ref. || (scalars and leptons only). 



results in three localized left- (right-) handed in four-dimensional sence zero modes of 
L (E) which describe three generations of usual leptons. The zero modes have the 
form 
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where r, 9 are polar coordinates in the bulk, p = 0, 1, 2 enumerates three zero modes 
which in the effective four-dimensional Lagrangian are described by two-component 
Weyl spinors \ p , e p , left- and right-handed, respectively. The radial functions L p (r), 
e p (r) have the following leading behaviour: 

l p (r) ~ r p , r -> 0; l p (r) ~ e~ 9lVr } r -> oo, 

and the similar for e p (r). 
The scalar potential 
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results in localisation of if and X inside the vortex. Interaction with the Higgs field, 



YHXL- — —E + YeH&L- — —E + h.c, 
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Figure 1: Sketch of the wave function profiles inside the brane. 



where Y and Ye are two Yukawa couplings, provides a hierarchical structure of masses 
and mixings of charged leptons, as is discussed in Refs. @, || . The reason for hierarchy 
is the different behaviour of three zero modes inside the brane which, after integrating 
out extra dimensions, results in the hierarchical structure of overlaps of wave functions 
for relatively narrow Higgs field (see the sketch at Fig. 1). 



3 Neutrino masses 

At the next step we follow the usual approach || to neutrino masses in models with 
extra dimensions, namely, we put a singlet fermion in the bulk. It plays the role 
of the usual right-handed neutrino of the see-saw models. The Dirac masses of the 
localized neutrinos are suppressed by the bulk size. The usual approach is to put 
three singlets in the bulk, otherwise only one linear combination of localized modes 
gets mass. This is in some contradiction with our intent to have a single generation 
in the higher dimensions. However, with more than one extra dimensions the bulk 
spinor has, in principle, enough components to give masses to more than one four- 
dimensional neutrino in the effective theory |12 |. Here, we exploit this fact in a slight 
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modification of our model of Ref. 0. 

Let us introduce the six- dimensional Dirac spinor N, which is a Standard Model 
singlet and has no half- integer axial charge under U(l) g . In this case, no interaction 
term similar to Eq. ([!]), which could result in localisation of the zero modes of N, can 
be written. The free spinor N satisfies the six-dimensional Dirac equation, 

iV A d A N = 0. 

It is convenient to perform the four-dimensional Fourier transform first, 

N = J d 4 ke- ik " x »N(K^r,9), 

substitute to the Dirac equation and multiply the latter by To, 

(k - kiToTi)N = DN, 

where the operator D = iToT a d a anticommutes with the operator in the left hand 
side. The latter fact allows us to look for a solution expanded in series of eigenvectors 
of D with eigenvalues uj, 

DN = tuN. 

The solution to this equation can be found in Appendix We note here that the 
zero mode uo = which is regular both in the origin and at infinity consists of just 
four linearly independent constant two-component spinors, 

/ ni \ 
n 2 

n 3 
V n 4 y 



N 



(3) 



From the four- dimensional point of view, we have two left-handed Weyl spinors n 2) 3 
and two right-handed Weyl spinors 11^4 (see discussion of chirality in six and four 
dimensions in ref. 0). The latter two independent spinors, indeed, can be used to 
give masses to two of three localized neutrinos. This provides three Am 2 which appear 
in the usual four- dimensional neutrino models. 

The advantage of having an extra Weyl spinor, however, disappears if we use 
interactions with H, X, and similar to Eq.(^). Indeed, U(l) g invariance requires 
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fields 


profiles 


charges 


representations 






u(X) 9 


U(1) Y 


SU{2) W 


SU{3) C 


fermion N 


bulk modes 


+1 





1 


1 


scalar S 


S(r) 

S(Q) = v s , S(oo) = 


+3/2 





1 


1 



Table 2: Fields introduced to give masses to the neutrinos. 

that N has the charge 3/2 under this group (half- integer axial charges are not allowed 
in order not to trap the iV particle on the brane) for the term like 

HXL - ~ Tj N 
2 

(or —3/2 if we use (1 + T- 7 )/2). Without the projector (1 ± r 7 )/2, the term is not 
gauge invariant; on the other hand, this projector kills one of two spinors n 14 . To 
overcome this difficulty, we introduce a new scalar field S which is analogous to X 
but has a charge 1/2 under U{l) g (the two fields, N and S, which we introduce in 
this paper, are described in Table 2Q). With these charge assignements, we can write 
down two interaction terms, 

yi HSL^±N, (4) 

and 

y 2 HS*L l -±^N. (5) 
Substituting zero modes of I and N, one gets the effective Lagrangian 

J rdrdO (y l HS(r)l 2 ^ p (r)e- i ^ 6 l p n 1 + y 2 H S (r)l p (r)e vp \n A ) 

Since H and S have no 9 dependence in their classical profiles, the integral over 9 is 

non-zero for p = 2 (the first term) and p = (the second). Thus, the modes 1^2 and 

^With these charge assignments, Dirac mass term for N is not forbidden. The simplest way to 
forbid it is to assign a mixed, axial and vector, charge to N. The axial charge qA should not be 
half-integer, otherwise N is localised on the brane. Together with the charge qs of S, they satisfy 
qs = 3/2 + qA, while the vector charge of N, qv = 1. For example, one of the solutions is qA = 1/4, 
qv = 1, qs = 7/4. 
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n4,i are paired in the mass matrix. The masses are given by integrals 

yi i2 27r JrdrHS(r)l (r)\n 4tl \ 

where the norm \m\ ~ 1/-R, R being the size of the extra dimensions. Masses of 1q 
and I2 differ only by Yukawa couplings yip in two terms. 

Two notes are in order. First, would we add another scalar with the same quantum 
numbers as S but with winding number 1, we would receive a more complicated mass 
matrix (mixing would appear in the neutrino sector as well as in the sector of charged 
leptons), but not the mass for the third neutrino - because we have only two right- 
handed Weyl spinors, ni and 114, available. One linear combination of l p would still 
be massless. Second, the Eqs. (f|), (||) mix \ p with nonzero modes of iV as well. 
However, to satisfy simultaneously constraints from weak charge universality and 
from Supernova 1987a, these mixings should be negligible. 

4 Parameters and estimates 

In order to count parameters and to estimate the scales involved in the model, we 
make use of the "step Higgs" approximation (full numerical analysis will be presented 
elsewhere fl3|]). First of all, recall that in six dimensions, a scalar field has dimension 
(mass) 2 and a fermion field has dimension (mass) 5 / 2 . This means that the trans- 
verse part of the fermion, described by l p (r), has dimension (mass) 1 (four-dimensional 
spinors, denoted by letters in bold face, are of course (mass) 3 / 2 ). Dimensions of cou- 
pling constants are [Y] = (mass) -3 , [g] = (mass) -5 , [h] = (mass) , [A] = (mass) -2 . 
The fundamental dimensionful parameters of the theory are the vortex scale v (in 
theories where large extra dimensions are invoked to reformulate the gauge hierarchy 
problem, v is of order mult i- dimensional Planck mass), and the compactification scale 
R. It is convenient to rescale the parameters of the theory: 

9i,e = gi, e v 5/2 } v = v 1/2 , p, = fiv~ 1/2 , v x = t>it>~ 1/2 , h = hv 1/2 , Y = Yv 3/2 , A = Xv. 

Tilded coupling constants g, h, Y, A are dimensionless, while v, vi, p, have dimension 
(mass) 1 . 



9 



In the narrow Higgs approximation, H(r) can be substituted by 

H(r) ~ 



H(0), r<r H - 
0, r > r H , 



where r# ~ (MO 1 is the width of if(r) which is supposed to be smaller than widths 
of other profiles. Then, we can use the following approximations for r < r#: 

X(r)«X(0), S(r)*S(p), 



The coefficients at r p in fermionic radial functions are determined, approximately, in 
Ref. fL3| . For the scalar fields if, X, and 5, the boundary conditions at r = are 
dH(r)/dr = 0, etc.; so that H(0), X(0), 5(0) are to be determined by solving the 
full nonlinear system [T3|]. We denote H(0) = X(0) = xvh, 5(0) = svh, where x 



and s are coefficients of order one, vh depends non-trivially on \i and v. Since /x and 
i> are independent parameters of the lagrangian, vh does not need to be very close to 
v, though we do not wish them to differ too much in order to avoid fine tuning of the 
parameters. 

Diagonal elements of the mass matrix of charged leptons are, to this approxima- 
tion, 



m pp « 2nYH(0)X(0) / r dr L(r)eJr) « vrFa;^ 

Jo u 



where 5 = \Jg{g e j h 2 should be sufficiently small to provide the mass hierarchy between 
three generations. Smallness of 5 exactly corresponds to the fact that H(r) is more 
narrow than l(r), e(r). 

Non-diagonal mass matrix elements are estimated in a similar way, taking into 
account the behaviour of F(r) at small r: 

vr 

The result is that 
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where generation-independent constant 




This corresponds to the mass matrix of charged leptons 

/ 1 \ 

M e ^M 0e i 5/2 

v lb 1 2 5 2 /3 j 



(6) 



where the overall constant Mq c = nxYvfj/v. 

It is straightforward to obtain, to the same approximation, the neutrino mass 
matrix 

/l 






M„ = Mi 



\ 








where 



M 



h 2 v vR 



and we rescaled, as before, 3/1,2 = yi^v 3 ^ 2 . The matrix is 3 x 2 because we have only 
two light right-handed modes, ni and 114. 

Thus, up to dimensionless constants of order one, the neutrino mass scale is sup- 
pressed by a factor of — with respect to the mass scale of charged fermions, v%/v. 

vR 

One can have simultaneously v 2 H /v ~ 100 GeV, m„ ~ 0.1 eV, and Planck mass 
relation v 2 R ~ 10 19 GeV satisfied for v ~ 1000 TeV, v H ~ 10 TeV, R ~ 10~ 6 mm. 

In our model, neutrino mass matrix does not contain mixings. In the sector of 
charged leptons, however, mixing appears, see Eq. (Q). Three parameters M 0e , 5, and 
e, determine the masses of three charged leptons; one can fit known values of M CifljT 
with M 0e ~ 650 MeV, 5 ~ 0.13, e ~ 2.5. The latter value means that mixing is not 
negligible in the leptonic sector. Since all leptonic mixings appear in Eq. (^j), our 
model predicts values of the neutrino mixing angles which are implicitly determined 
by masses of charged leptons. In usual notations, they correspond to sin 8 12 ~ 0.14, 
sin# 13 ~ 0.37, sin ^23 ~ 0.99. On the other hand, values of two Am 2 are predicted 
only by the order of magnitude since they depend on unknown constants yi^. 
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5 Conclusions 



In a class of multidimensional models with one vector-like fermionic family, the low- 
energy effective theory describes three chiral families in four dimensions. Hierarchy 
of fermionic masses appears without fine tuning of parameters. In this paper, we 
have shown that nondegenerate neutrino masses can be incorporated in a model of 
this class with the help of only one multidimensional Dirac fermion. One of three 
neutrinos is exactly massless| while masses of two others are parametrically close to 
each other. Neutrino mass matrix is diagonal, but mixings in the leptonic sector are 
present in the mass matrix of charged leptons. Once masses of e, fi, r are fixed to 
their experimental values, the model predicts neutrino mixing angles, sin #12 ~ 0.54, 
sin# 13 ~ 0.31, sin 6*23 ~ 0.99. Of course, these values were obtained with very rough 
approximations, and should not be taken too seriously. Neutrino mass differences 
are predicted by order of magnitude only, Am 2 ~ (0.1 eV) 2 for natural values of 
parameters. 
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I.Gogoladze, A.Neronov, E.Nugaev, V.Rubakov, G.Senjanovic, M.Shaposhnikov, 
S.Sibiryakov, A.Smirnov for helpful discussions on different subjects related to this 
work. The work of M.L. and S.T. is supported in part by RFFI grant 99-02-18410a, 
by the Council of Presidential Grants and State Support of Leading Scientific Schools, 
grant 00-15-96626, by CRDF award RP1-2103, and by the programme SCOPES of 
the Swiss National Science Foundation, project No. 7SUPJ062239, financed by Fed- 
eral Department of Foreign affairs. The work of S.T. is supported in part by Swiss 
Science Foundation, grant 21-58947.99. This work was initiated during the visits of 
M.L. and S.T. to ULB, which we thank for kind hospitality and partial support by 
the 'Actions de Recherche Concretees" of "Communaute Frangaise de Belgique" and 
IISN-Belgium. S.T. thanks High Energy Group of ICTP (Trieste) for hospitality dur- 
ing his visit, when part of this work was done and numerous discussions were carried 
out. 



§See, however, Appendix 
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A Solution to the Dirac equation for N 



In this Appendix, we present the solution to the transverse part of Dirac equation, 



DN = wiV, 



where D = iT T a d a . The general solution for to = is 



N 



-pi 



JpiO \ 



A/2 2 e 



Here, n a are independent two-component spinors. The solutions regular both at 
small and large r correspond to p a = 0, Eq. (||]). For them, the normalization factors 
Na — V^R, where R is the radius of the bulk. 

For u ^ 0, the regular at r = solution depends on two two-component spinors 
n, m for each u: 

( n UJ J p (tor)e ipe \ 
m u) J q (ur)e iqe 

A series of eigenvalues 00^ are determined by the boundary conditions at large r. 
The exact eigenvalues depend on the compactification scheme chosen (see discussion 
in Appendix 0), but in general, they are of order 



N 



(7) 



to, 



(p) 



nn/ R. 



B Compactification 

The compactification scheme defines boundary conditions at large r. It should satisfy 
the following criteria: 

(i) it should preserve three zero modes with different winding numbers of the fields 
localized on the vortex; 

(ii) it should allow for at least one (Dirac) zero mode of the bulk singlet. 
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The second requirement is essentially equivalent to compactification on a Ricci-flat 
manifold (that is, a manifold with zero scalar curvature TV) fl4j . Toroidal compact- 



ifications satisfy this requirement since the torus have 1Z = 0. However, the corre- 
sponding boundary conditions usually kill one or two of three zero modes because 
they allow either even or odd winding numbers for a given field. Compactification of 
a vortex on a sphere preserves three zero modes |15j but does not satisfy the second 
requirement: for a sphere, TZ ^ 0. This would open an interesting possibility of hav- 
ing three massless neutrinos and explaining observed neutrino anomalies by mixing 
with Kaluza-Klein modes only (this would lead to a nonconventional pattern of os- 
cillations, as was discussed in Ref. |16[ in the frameworks of a toy model). However, 
this possibility is in contradiction with bounds on mixing with Kaluza-Klein states 
from supernova evolution. 

One can also keep the compactification mechanism unspecified, and consider flat 
space of finite volume. This means that the fields will live inside a cilinder of large 
radius R in six dimensions. The corresponding boundary condition, that is, absence 
of fermionic current through the boundary, can be reformulated WR as 



iT A d A ip = V> 

for a Dirac spinor ip. Applied to N, Eq. ([7]), this boundary condition selects not only 
constant, but also rather peculiar zero modes which grow up at large r. In this case, it 
is not necessary to introduce the additional field S to obtain a second neutrino mass, 
and it suffices to use the additional mode, alternatively, with the field S included, 
3 (even 4 if the model is extended) neutrino masses can be generated, although the 
last 2 are then hierarchically smaller than the ones previously defined. At this level 
however, such modes appear more as a quirk of a particular way to implement the 
finite volume. 

It seems the most promising to have a finite volume in the case of warped trans- 
verse dimensions. Particularly interesting is the fact that there exist six-dimensional 



solutions with gravity localized on an Abelian vortex, Ref. These six-dimensional 
solutions provide also a natural framework for the localization of gauge fields at weak 
coupling, thus opening the way to fully realistic and calculable models. 



14 



References 



[1] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D61 (2000) 033005 |Kep- 
ph/99034l7| ; 



G. Dvali and M. Shifman, Phys. Lett. B475 (2000) 295 ||hep-ph/000l"072 
T. Gherghetta and A. Pomarol, |hep-ph/ 0003129] 



D. E. Kaplan and T. M. Tait, JHEP 0006 (2000) 020 [ |hep-ph/00042U0 
S. J. Huber and Q. Shan, frep-ph/0010195l 



[2] M. V. Libanov and S. V. Troitsky, Nucl. Phys. B 599, 319 (2001) fKep- 
ph/0011095 



[3] J. M. Frere, M. V. Libanov and S. V. Troitsky, Phys. Lett. B 512, 169 (2001) 
|frep-ph/0012306 . 



[4] E. J. Weinberg, Phys. Rev. D24 (1981) 2669; 

[5] N. Arkani-Hamed, S. Dimopoulos, G. Dvali and J. March-Russell, |hep- 



ph/9811448 



K. R. Dienes, E. Dudas and T. Gherghetta, Nucl. Phys. B 557 (1999) 25 jhcp- 



ph/9811428 



[6] G. Barenboim, G. C. Branco, A. de Gouvea and M. N. Rebelo, hcp-ph/0104312 . 
[7] A. De Gouvea, G. F. Giudice, A. Strumia and K. Tobe, [hep-ph/0107156 . 



[8] R. Barbieri, P. Creminelli and A. Strumia, Nucl. Phys. B 585 (2000) 28 jhep- 
ph/ 0002199| ; 



A. Lukas, P. Ramond, A. Romanino and G. G. Ross, JHEP 0104 (2001) 010 
|hep-ph/0011295 ]. 

[9] G. Dvali and M. Shifman, Phys. Lett. B 396 (1997) 64 [Erratum-ibid. B 407 
(1997) 452] [|hep-th/9612128jl . 



[10] E. K. Akhmedov, |hep-th/0 107223 



[11] S. L. Dubovsky and V. A. Rubakov, |hep-th/ 01 05243 



15 



[12] K. Agashe and G. Wu, Phys. Lett. B 498 (2001) 230 [|hep-ph/0010117 



[13] E.Ya.Nugaev, Master Thesis from the Moscow Institute of Physics and Technol- 
ogy, 2001; M.V.Libanov and E.Ya.Nugaev, to appear. 

[14] E. Witten, "Fermion Quantum Numbers In Kaluza-Klein Theory," PRINT-83- 
1056 (Princeton). In: Appelquist, T. (Ed.) etai: Modern Kaluza-Klein Theories, 
438-511. 

[15] A. Nakamura, S. Hirenzaki and K. Shiraishi, Nucl. Phys. B 339 (1990) 533; 
A. Nakamura and K. Shiraishi, Mod. Phys. Lett. A 5 (1990) 1109. 

[16] N. Cosme, J. M. Frere, Y. Gouverneur, F. S. Ling, D. Monderen and V. Van 
Elewyck, Phys. Rev. D 63 (2001) 113018 ||hep-ph/ 00 10192^ . 

[17] K. Johnson, Acta Phys. Polon. B 6 (1975) 865. 



[18] T. Gherghetta and M. Shaposhnikov, Phys. Rev. Lett. 85 (2000) 240 jhc^ 
| th/ 00040141 ; 

M. Giovannini, H. Meyer and M. Shaposhnikov, [hep-th/0104118 . 



16 



